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Abstract 

We discuss a model for a universe with discrete matter content instead of the continuous perfect fluid 
taken in FRW models. We show how the redshift in such a universe deviates from the corresponding 
one in an FRW cosmology. This illustrates the fact that averaging the matter content in a universe 
and then evolving it in time, is not the same as evolving a universe with discrete matter content. The 
main reason for such deviation is the fact that the photons in such a universe mainly travel in an 
empty space rather than the continuous perfect fluid in FRW geometry. 
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Chapter 1 



Introduction 



The main goal for this work is to show a model for the universe in which it is not assumed that 
the matter content of the universe is, on average, a perfect continuous fluid which permeates all the 
universe. We discuss a more realistic model in which the matter content is assumed to be discrete and 
show how this does not affect the dynamics of the universe but affects the light propagation in the 
universe. 

The main motivation behind discussing such model is that in FRW cosmology, the most usable 
model for the universe, the matter content is assumed to be a perfect continuous fluid which permeates 
the universe. To use this model to interpret the data, observed experimentally, we should assume the 
existence of an exotic continuous perfect fluid with a constant density through all the universe i.e. it 
never dilutes as the universe expand. This perfect fluid is called dark energy. Therefore, considering 
such an inhomogeneous model could provide further insight to investigate the validity of assuming the 
existence of dark energy. 

We start off in Chapter 2 by discussing Einstein's picture of gravity and the geometric consequences 
for this approach; we discuss the geometric description of the Einstein field equations in a general 
curved spacetime. We then discuss the general form of the energy-momentum tensor and then write 
it explicitly for the perfect fluid case in a general frame of reference. 

In Chapter 3, we discuss two exact solutions for Einstein's field equations, namely, the Schwarzschild 
solution, and the Friedmann-Robertson- Walker (FRW) solution. In each geometry, we discuss the 
equations of motion of photons and massive particles in the resulting spacetime and also the redshift. 1 
In the FRW solution, we get two independent equations for the scale factor of the universe. We then 
discuss the three solutions in the case of the matter-dominated universe, namely open, closed and flat 
universes in that case. 

In Chapter 4, we use this quick introduction to discuss a discrete model for the universe introduced 
by Lindquist and Wheeler (LW) in 1957 in a seminal paper [9]. They discussed in their model the 
case of a matter-dominated closed universe and showed that the global dynamics are similar to the 
corresponding case in the FRW cosmology. We then discuss how this model was extended in another 
wonderful paper by T. Clifton and P. Ferreira in 2009 [2] to discuss the light propagation for such a 
set-up. In the case of a flat universe, they showed a deviation in the redshift from the corresponding 
redshift in an FRW cosmology, namely, the Einstein-de-Sitter (EdS) model. The main advantage of 
such a model is that that global dynamics of the FRW cosmology was shown to be a limiting case of 
that model, namely, when the number of the discrete islands becomes very large. 

It is also important for the reader to notice that in the whole of this work, we are taking the metric 
with a —2 signature. Otherwise stated, the sign of the time part is positive and the sign of the spatial 
part is negative, i.e. in flat spacetime the metric is expressed as ds 2 = dt 2 — dx 2 — dy 2 — dz 2 . 



lr The redshift is one of the most important quantities in observational cosmology. 
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Chapter 2 

Einstein Theory for Gravity 



2.1 A different approach to treat the gravitational force 

Two important observations were a cornerstone, and a main motivation for Einstein in developing the 
General Theory of Relativity, a theory for gravity. First, the experimental fact that the inertial mass 1 
of a massive particle is the same as its gravitational mass. 2 And second, the Equivalence Principle 
which says that in a freely falling and non-rotating part of space, i.e. an elevator, the laws of physics are 
the same as the special relativity laws. Otherwise stated, if we consider only a uniform gravitational 
field in a closed elevator, i.e we ignore the tidal forces, we can't distinguish between the gravitational 
field and the field which would be developed due to a uniform acceleration for that part of space. 
This means that all the particles in that region of space have no acceleration in the inertial frame of 
reference moving with this region of space. 

This led Einstein to think of gravity as a curvature in the spacetime which is caused by matter and 
not as a force. Einstein's main idea was to replace the equation of motion for a particle moving under 
gravitational forces by the equation of motion for a freely-moving particle in a curved space in which 
the curvature depends on the amount of matter and energy in it, i.e. 

Here, P is the four-momentum, and t is the proper time. 

This means that the spacetime trajectory of the freely-moving particle in curved spacetime is a 
geodesic of that spacetime. 

2.2 Einstein field equations of gravity 

They are a set of 10- independent, non-linear differential equations for the metric fields g fll/ (x fJ '), 3 they 
relate the spacetime curvature to all the sources of energy in it, i.e. ordinary matter, radiation, and 
dark energy. They can be written in tensorial form as follows. 

• With a zero cosmological constant, A, i.e. the amount of the dark energy is zero, 

R^v - ^Rg^u = —nTpv. (2.1) 

• With a non-zero cosmological constant, A, 

Rfj.v — ^Rg^u — —kT^u — Ag M „. (2.2) 

Here, R uv is the Ricci tensor, g M „ is the metric tensor, R is the curvature scalar, and n = — =- with 
the Newtonian gravitational constant G. 



1 Inertial mass represents the ability of a massive particle to resist any change in its state of motion. 

2 Gravitational mass represents the ability of a massive particle to accelerate toward other fixed gravitating object. 

3 Throughout this essay we consider /i, v, <5, A, p, a, a, b, c,d,e = 0, 1, 2, 3 (spacetime), but i, j, k, I = 1, 2, 3 (space only). 



4 



It is common to call the right-hand side of Einstein held equations as the Einstein tensor = 
Also, R a b and g a b are related by 

-Rah = R C abc = ^°kb^ ac — ^°kc^ ab + dbT c ac — d c T c a b, (2.3) 

where R d a bc is the general curvature tensor, or the Riemann-Christoffel tensor, and T b ac is the afhne 
connection. 

The afhne connection is related to the metric elements by 

r\ c = \g ad {d b g dc + d c g db - d d g bc ). (2.4) 

Equations (2.3) and (2.4) say that the Ricci tensor components depend on the metric components g^ 
and their derivatives. Therefore, we measure the spacetime curvature 4 using g^ v . This is the reason 
for considering the held equations as a differential equations on the metric fields g^v 
If we contract equations (2.1) and (2.2) with g^, then, 

• from equation (2.1), we get R = kT, therefore the Einstein held equations, with zero cosmological 
constant, can be rewritten as 

r 



R^v = -k(T^ v - -Tg^ u ). (2.5) 



• from equation (2.2), we get R = kT + 4A, therefore the Einstein field equations, with non-zero 
cosmological constant, can be rewritten as 

Rftu = -iztTfja, - -Tg^u) + Ag^. (2.6) 

Here, T = g^T^" = Tfc. 

Now we can write equations equivalent to the equations of motion of particles moving under the 
influence of a gravitational field. These will be the equations of the geodesic 5 in a general curved 
spacetime, if that spacetime is curved due to the sources of that gravitational held. If the general 
trajectory in the spacetime is cc' i (s), where s is the parameterization parameter along that curve, the 
tangent for that curve will be 6 

dx^ 

u » = x » = ——, 
as 

and for that curve to be a geodesic, the covariant derivative of x^ should vanish, so we write the 
geodesic equations as 7 

V = + V u° = ^ + = 0. (2.7) 

The full trajectory for the particle, i.e. the solution of the equations of motion, is then given by the 
functions x M (s), that obey equation (2.7). 

Provided that we accept the new approach for gravity, we can conclude from the above discussion 
that the metric components of that curved spacetime play an essential role in determining the trajecto- 
ries of particles moving under the effect of gravity. They depend only on the energies in the spacetime 
which cause the spacetime to curve through the Einstein held equations. 



2.3 Description of the Einstein field equations 

2.3.1 The spacetime curvature-sources side of the equation 

The right-hand side of the equation (2.2) contains two different kinds of energy: the energy due to 
matter and radiation, in T^ v ] and the cosmological constant, which expresses the dark/vacuum energy, 
i.e. the energy in the universe if we remove all the matter and the radiation from it. 8 

4 The spacetime deviation from being flat. 

5 Geodesic line is the shortest possible line between two points in a general curved spacetime. 

6 Here, the dot indicates the derivative with respect to s. 

7 For a more detailed mathematical description the reader is referred to [7]. 

8 In other words, one can think of the dark energy as the ground state of a quantized-bounded system, which is 
normally not of zero energy because of the Heisenberg uncertainty principle. 
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The energy- momentum tensor, T Ml/ , 9 is the term which expresses the energy and the momentum 
flow for radiation and matter, as a source of the gravitational field that curves the spacetime, in general 
it can be expressed as 

T^(x^) = pix^u". (2.8) 

Here, u M = u^(x^) = j u (c, u(x^)) and p is the energy density. 

All the above quantities, p, T Miy , and u M are measured in the frame of reference described by the 
coordinate system x^ and at a certain point in that spacetime. Therefore, at that point 

• T 00 = expresses the energy density of the gravitating particle, 

• T at — j^pcu 1 relates to the energy flux in the i th direction of the space, 

• T l ° — ry^pcu 1 relates to the momentum flux in the i th direction of the space, and 

• T zj = juptfuJ represents the rate of the flow of the i th component of the momentum in the j th 
direction. 

Conservation of the energy and momentum can be expressed in the general spacetime as 

V^T"" = d^ v + T» x T Xu = 0, (2.9) 
where V M is the covariant derivative operator. 



2.3.2 The Spacetime curvature side of the equation 

The left-hand side of equation (2.1), or (2.2), measures how much the spacetime is curved or deformed 
from being flat due to the existence of the deformation sources, i.e. any type of energy mentioned 
above. The Ricci tensor measures how much the spacetime, a 4-dimensional manifold, deviates from 
being a 4-dimensional Euclidean space. As indicated in the previous section, the determining quantities 
on this side of the equation are the metric tensor quantities, g^ v . 

In the literature, there are many routes to derive the Einstein field equations. One of them, as 
mentioned before, is to use the conservation of energy condition (2.9) and the fact that these equations 
should, in non-relativistic limit, be reduced to the Newtonian field equations, to restrict the possible 
forms of the equations, for more details about the derivation the reader is referred to [7], or to [1]. 

Another way to derive them is to take the Einstein-Hilbert action S = Ry/\g\d 4 x, and apply the 
variational principle, the principle of least action, to get the field equations. For more mathematical 
details, the reader is referred to [4] section 37.4. 



2.4 The energy momentum tensor in some special cases 
2.4.1 For a free space 

The energy-momentum tensor components for that part of spacetime vanish since it contains no matter, 
i.e. T MI/ = T = 0. Assuming no dark energy, this gives, from equation (2.5), that i? M „ = 0. In presence 
of dark energy, this gives, from equation (2.6), that — Ag^. 



2.4.2 For a perfect fluid 



The advantage of taking the perfect fluid case is that in this approximation there is no interaction 
between the particles, i.e. the sheer forces are ignored, and the thermal motion of the particles can 
be ignored. Using these assumptions, we can write T for a perfect fluid in its instantaneous frame of 
reference (IFR) as 



T = [T^] 



' 1 

c p 
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p 














p 



(2.10) 



^T^ v is called sometimes the stress-energy tensor 



G 



where p is the pressure of the fluid, and p is the energy density of the fluid. It can also be written in 
tensorial form as 

T" u = {p + p/c 2 )u»u" -prf v , 
where rf" is the Minkowski metric which is given by 

1 ' 
0-10 

0-10 

-1 



n 



{IV 



Therefore, the full tensorial expression for T in any arbitrary coordinate system is 

T^ v = (p+p/c 2 )uV -pgT- 



(2.11) 



If we substitute this T M ", given in equation (2.11), in the covariant expression of the conservation of 
energy, given by equation (2.9), we get the relativistic continuity equation and the relativistic equation 
of motion of a perfect fluid. Namely, in flat spacetime it becomes d^T^ v = 0, which gives, using 



u^Ufj, = c 2 , 



• the continuity equation 



d^ifm") + (p/c 2 )d^ = 0; 



• and the equations of motion 

( P + P /c 2 )(d^)u^ = (rr - u^/c 2 )d^p. 



(2.12) 



(2.13) 
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Chapter 3 

Exact Solutions for Einstein Field 
Equations 

3.1 Introduction 

Accepting Einstein's new approach for treating the gravitational forces, we present here two exact 
solutions for his field equations. Namely, the Schwarzschild solution, 1 which is the spacetime geometry 
outside a spherically symmetric matter distribution, and the Friedmann-Robertson- Walker (FRW) 
solution, which is the spacetime geometry of a homogeneous and isotropic space containing a uniform 
distribution of matter. 

3.2 Schwarzschild geometry 

Since the space outside that spherically symmetric mass distribution is free, 2 therefore, as we mentioned 
in Section 2.4, 

= 0. (3.1) 
The general form of a static, spherically symmetric metric in polar-spatial coordinates is 

ds 2 = A(r)dt 2 - B(r)dr 2 - r 2 d9 2 - r 2 sin 2 Odcj) 2 . (3.2) 

Using condition (3.1) and the line element (3.2) we can solve for A(r), and B(r) and we then get the 
Schwarzschild line element, 3 written as 

dr 2 

ds 2 = c 2 (l - 2m/r)dt 2 - — - r 2 d9 2 - r 2 sin 2 6d<j) 2 . (3.3) 

(1 — 2m/r) 

Here, m = MG/c 2 , where M is the total mass of the spherically symmetric mass distribution, r is the 
Schwarzschild radial component, 4 t is the Schwarzschild time, and and <j> are the polar angles. 

By looking to the second term in equation (3.3), we can see that r is not the same as the radial 
distance in normal spherical polar coordinates. Also, from this expression for the metric we can 
interpret t as the proper time experienced by a stationary observer at r — > oo. 

It was shown later by Birkhoff that the spherically symmetric vacuum solution is the same for a 
static and non-static spacetime 5 , since both have the same line element (Birkhoff theorem). 

In the vicinity of this mass distribution M, following Einstein's proposal, the spacetime will be 
curved, and to determine how other objects, photons or massive particles, behave, we should determine 
the geodesic equations for a freely-moving object in a spacetime with the line element given by equation 
(3.3). 

1 Schwarzschild was in the trenches on the Eastern Front conflict when he came up with his solution, and sadly he 
did not survive the conflict. 

2 It is important here for the reader to notice that we shall ignore the dark energy term in our treatment. 
3 For a detailed derivation the reader is referred to [8] pages 142-7. 

4 r = 2m is dubbed the Schwarzschild radius, and if the massive object contract within r = 2m, this object becomes 
Schwarzschild blackhole. 

5 Here, we taking about the spacetime exterior to the matter distribution. 
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3.2.1 The red-shift in the Schwarzschild geometry 

For the case where both the emitter (E), and the receiver (R) are stationary in the space, 6 the red-shift 
z is expressed as 

' 1 — 2m/ Vr 




l-2m/r,' (3 - 4) 

where te and tr are the Schwarzschild radial coordinates for the emitter and the receiver, respectively. 
For a detailed proof of equation (3.4), the reader is referred to read the elegant proof in section 4.3 of 
[5]. 

3.2.2 The geodesic equations in the Schwarzschild geometry 

We derive the geodesic equations using the Euler-Lagrange method. ' The Lagrangian L is given by 

f 2 

L = g^x" = c 2 (l - 2ni/r)i 2 - - — — ^ - r 2 (6 2 + sin 2 6»</> 2 ). 

i , ^ i T ■ d ( dL\ dL 

if we apply the Euler-Lagrange equation, — — — = — — , we get 

as \ox^ J oxf 1 

(1 - 2m/r)i = k, (3.5) 

"'"V- 2 . ™ . . 2 r 2 -r(0 2 + sin 2 ^ 2 )=O, (3.6) 



1 — 2m/ V r 2 r 2 (l — 2m/ r) 2 
■a 2 . 



r# — sin cos 6>ti 2 



r 



(3.7) 



r 2 sin 2 00 = h, (3.8) 

where fe, and h are constants. The moving particle total energy is E = in particular, for a 

particle at rest at infinite, i.e. f = 0, and r — oo, its 4-velocity will be u M = (1, 0, 0, 0). Therefore, 

2 ^ 
E = po — goQiriot = kmoc =^> fc = 



TO C Z 



where mo is the rest mass of the moving particle. Directly from equation (3.8) we see that h = — I, 
where I is the angular momentum I of the particle. 8 

We then consider the particles moving in an equatorial plane 9 , i.e. we keep only solutions with 
9 = 7r/2. This satisfies equation (3.7), and simplifies the equations into 

(1 - 2m/r)i = k, (3.9) 



i lie 2 - 2 m 



i .., 2 r~ 2n 9 , ^ 2 r 2 - = 0, (3.10) 

r 2 $ = h. (3.11) 
Also, we can simplify these equations further by considering one of the following constraints: 

• g^ii'X^x 1 ' = c 2 for massive particles. 

• gtiv&^i 1 ' = for photons. 

We then replace the second-order differential equation (3.10) by one of these constraints, a first-order 
differential equation. Now we have the geodesic equations, and by solving them, we get the functions 
x^(s) that represent the trajectory for the particle. Here we separate the two cases: massive particles, 
and photons. 



6 Stationary here means that r, 8 and <j> are of constant values. 
7 The following treatment is based on the treatment on [7]. 
s l = g 33 j> = -r 2 sin 2 d<t> = -h. 

9 This idea of taking the equatorial plane will only simplify the calculations, with no loss of generality since it is always 
possible to rotate the spatial part of the Schwarzschild coordinates such that the trajectory of the particle lies in the 
(r, <f>) plane. 
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1. The trajectories for massive particles. 

Replacing equation (3.10) by the condition g (lv x^x l/ — c 2 , we get 

(l-2m/r)i = k, (3.12) 
p. 

c 2 (l~2m/r)i 2 — --r 2 2 =c 2 , (3.13) 

1 — 2m /r 

r 2 <P = h. (3.14) 
By substituting equations (3.14) and (3.12) into equation (3.13), we get the energy equation 

f 2 + ^(l-2m/r)-^=c 2 (fc 2 -l). (3.15) 
dv h dv 

Using equation (3.15) and the fact that — = —^-rj, with u = 1/r, we get the it-equation, which 

dr r z d(f> 

determines the shape of the massive particles' orbits, 

d 2 u GM 3GM 2 . , 

d^ +u =J^ + — «■ (3-16) 

To investigate these two equations (3.16) and (3.15) more, it is useful to consider two cases: 
radial motion and circular motion of the massive particles. 

• Massive particles moving radially. 

We have <j> = constant =>0 = O=>/i = O, and hence, using MG — mc 2 , the energy equation 
becomes 

•2 2/, 2 ^ 2 MG .. MG 
r 2 = c 2 (k 2 - 1) H — =>■ r = ~ r . 

For a particle dropped from rest at r = R c 2 (k 2 - 1) = - 2J j^ r 2 = 2MG (i - i). 
If the particle is dropped from rest at infinity, i.e. R = oo, we get that i = 1 _ 2 m/r an( i 5 

(a) for the in-falling particles f — — \j 2M r G ", and 

(b) for the out-falling particles f — : 1 ' ' ' 



For the in- falling case, we have = — J 2M r G , which by integration gives 



2 



3 V 2mc 2 



(3.17) 



where ro is the Schwarzschild radial coordinate at r = 0. To find the trajectories in the 
(r, <)-plane, we use = ^ = — (1 — 2m /r) y / '2mc 2 /r , which by integration gives 



2 / / 7q / r 3 \ 4m ( r$ r 



3 V V 2mc 2 V 2mc 2 / c \ V % m V 2™ / 



2m 
In 

c 



Vr/2m+ l \ ( y/r /2m-l 
v /r/2m - 1 J \y/r /2m + l 



We notice from the above two equations that at r = 0, r is finite, but, at r = 2m > 0, t 
goes to infinity, which means that the particle takes finite proper time to reach r = 0, but 
for a stationary observer at r — oo, it takes infinite time before reaching r — 2m. Figure 
3.1 shows the radially in-falling trajectory of a particle released from rest at infinity. 



1 This figure is taken from [7], and in the figure fj, is the same as our m. 
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Figure 3.1: Trajectory of a radially in-falling particle released from rest. The dots correspond to unit 
intervals of cr/m, where r is the particle's proper time and we have taken r = t = at = 8m. 11 



Massive particles moving in a circular orbit. 

d 2 u 

Here we have that r — constant, and so r = f = 0, also u = constant or —pr 
consequently, the u-equation gives 



K 2 = 



MGr 2 
r — 3m 



and the energy equation gives 



E 1 - 2m/r 



m c 2 y/\ - 3 TO / 7 



(3.18) 



(3.19) 



Equation (3.19) shows that the particles' orbits with radial coordinate r > 4m are bounded, 
i.e. k < 1, the orbits with 3m < r < 4m are not bounded, and for a particle to orbit at 
r = 3m it needs to have infinite energy. This is shown in Figure 3.2. 

2. The trajectories for photons. 

Replacing equation (3.10) by the condition g^x^x" — 0, we get 12 

(1 - 2m/r)i = k, 

V.2 



c 2 {l - 2m/r)t 2 



r 2 6 2 = 0, 



(3.20) 
(3.21) 



1 — 2m/r 

r 2 $ = h. (3.22) 
By substituting equations (3.22) and (3.20) into equation (3.21), we get the energy equation 



r 2 + -^(1 - 2m/ V) = c 2 k 2 . 



(3.23) 



12 For photons trajectories we can't use the proper time, r, as a parameterization parameter, so we use an arbitrary 
parameter s. 
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Figure 3.2: This plot shows the dependence of k = S total /moc 2 on the radial component for that orbit. 
It shows that for a particle to be orbiting at r = 3m, its total energy should be infinite, and the 
innermost stable orbit, i.e. with the lowest energy, is at r = 6m. 



dv h dv 

As for the massive particles case, using equation (3.15) and the fact that — = —s — , we can 

or r z dcp 

write the u-equation as 



d 2 u 



3GM 



(3.24) 



To investigate these two equations (3.24) and (3.23) more, it is useful to consider the two cases 
of radial motion and circular motion of the photons. 



• Photos moving radially. 

We have (j) — constant =>■ <f) = 



h = 0, hence, the energy equation becomes 
r = ±ck. 13 



(3.25) 



Equation (3.25) determines the trajectory in the (r, s)-plane, using it, with equation (3.20), 
we can write 

dv T 

— = t = ±c(l - 2m/r). (3.26) 

dc £ 

At r — > co, from equation (3.26), we can write that 4| = ±c, which represents the normal 
light cone with slope ±1. 

From equations (3.28), and (3.27), when r — > 2m =>■ the slope of the light cone — > ± oo, 
which means that the light cone closes up when the Schwarzschild radial coordinate ap- 
proaches 2m, and that explains the result that the massive particles take infinite Schwarzschild 
time to approach the Schwarzschild radius r = 2m. Figure 3.3 shows the structure of the 
lightcone in Schwarzschild's geometry, and how they close up when r approachs 2m. 
Using equation (3.26) we can learn about the trajectory in the (r, t)-plane, namely, 
— for in-falling photons, 



ct 



2m In \ r/2m — 1| + constant. 



(3.27) 



13 Here the + and the - signs correspond to the out-falling and the in-falling photos, respectively. 
14 This figure is taken from [7], and in the figure fi is the same as our m. 
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r=0 r = 2ft 



Figure 3.3: Lightcone structure of the Schwarzschild solution. 1 
— for out-falling photons, 

ct = +r + 2mln|r/2m- 1| + constant. (3.28) 
• Photons moving in a circular orbits. 

d 2 u 

Here we have that r =constant r = f = 0, and u = constant =>■ = 0. Consequently, 

cup 

the it-equation gives that the only possible circular photons' orbit is at r — 3M 2 . 

For the sun, 3A i ~ 4.5 km, but the sun radius is much larger than that, so we can't see 
such circular orbit for photons around the sun. However, in most cases, that is possible for 
black-holes. 

If we calculate the Riemann tensor components, Rijkh for the Schwarzschild line element (3.3), we find 
some non-zero components which depend only on the Schwarzschild radial coordinate r. This shows 
that the spacetime geometry around a spherical mass distribution M is not uniformly curved, but still 
spherically symmetric since these non-zero components don't depend on the polar angles. 



3.3 Friedmann-Robertson- Walker geometry 

The main reason for studying the geometry of isotropic and homogeneous space is that, on a large 
scale, the universe looks isotropic, which is also bolstered by the constancy of the cosmic microwave 
background radiation temperature. It is also believed that there is no preferred center for the universe, 
i.e. there is no spacial place in the universe. This implies that the universe is also homogeneous. 

3.3.1 The metric for an isotropic and homogeneous space 

To describe a metric, contains such properties, we slice the spacetime into a continuous series of 3D 
spacelike 1 5 hyper-surfaces . 1 6 

To define global time, we slice the spacetime such that the hyper-surfaces form a non-intersecting 
spacelike continuous series, i.e. a constant value of the parameter t is assigned for each hyper-surface. 
Normally we construct these series, which can be constructed in many ways, such that the worldline 
for a comoving 17 observer is perpendicular to the hyper-surface, which the observer belongs to. This 
is shown in Figure 3.4. 

15 In general, in a spacelike hyper— surface any 2 events are separated by a spacelike distance in the spacetime, so we 
can always find an incrtial frame of reference in which the 2 events occur at the same time, which is so called the surface 
of simultaneity. 

16 In general, hyper— surface refers to an n-d manifold embedded in n + 1 Euclidean space. 

17 Comoving usually means that as time evolves for that observer its spatial coordinate values, (r, 9, <f>), remain constant. 
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For such slicing of spacetime, the worldlines are orthogonal to hyper-surfaces, which represent the 
spatial part of the spacetime. We can write the line element as 18 

ds 2 = c 2 dt 2 - g ij (t,x 1 ,x 2 ,x 3 )dx i dx j . (3.29) 

It is also worth noting that here the parameterization parameter, t, of the hyper-surfaces is propor- 



S pacelike 
hypers urface 




Simultaneity surface 
in local Lorentz frame 



Figure 3.4: Hyper-surfaces of simultaneity, to which the worldlines are orthogonal 
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tional to the proper time measured by a comoving observer. And that the worldline for such space, 
given by x M = (t, A, B,C), where A, £?,and C are constants, satisfies the geodesic equations (2.7). 

The metric given by equation (3.29) should represent homogeneous and isotropic spacetime. There- 
fore, we can rewrite it as 

ds 2 = c 2 dt 2 - S 2 (t)h. lj dx l dx j . (3.30) 

Here, S(t) is a scale factor which depends only on t 7 and hij are functions of the spatial coordinates 
x 1 only. We can rewrite the spatial part, da 2 = hijdx l dx^ , in terms of spherical polar coordinates as 

da 2 = B(r)dr 2 + r 2 (d9 2 + sin 2 9d(f> 2 ). 



(3.31) 



The 3D hyper-surface should be maximally symmetric, therefore, the curvature of the space should 
depend only on one quantity, K, and for such maximally symmetric space, we can write the Ricci 
tensor as 

Rij = -2K gij . (3.32) 

And, if we require the metric (3.31) to satisfy the condition (3.32), we find that B(r) — 1/(1 — Kr 2 ), 
so our spacetime metric can be rewritten as 

' h ' 2 +r 2 (d0 2 + sin 2 6d<p 2 ) . 



ds 2 = c 2 dt 2 - S 2 (t) 

1 — Ar 

We then make the following change of variables r - 

dr 2 



(3.33) 



ds 2 



where R(t) is given by 



(fdt 1 - R 2 (t) 



R(t) 



1 - kr 2 



S(t) 



K 

K\ l ' 2 r, and k = t—;, which gives 
is 



r 2 (d6 2 



(3.34) 



if K ^ 0, 



|if|V2 
S(t) if if = 



18 For more details about the following derivations the reader is referred to 
19 This picture is taken from [7]. 
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Here, r is an arbitrary radial coordinate, and the coordinates r, 9 and <p are comoving coordinates, i.e. 
as the universe evolves in time, the worldline for any galaxy, ignoring any peculiar velocity, has fixed 
values of (r,0,<j>). This is because, as we said before, the worldlines satisfy the geodesic equations. 

In the metric (3.34), k is 1,0 or — 1, which corresponds to close, flat, or open universe, respectively. 
Investigating these cases separately, we can rewrite the metric (3.34) as 

ds 2 = c 2 dt 2 - R 2 (t)[d X 2 + S 2 (x)(d6 2 + sin 2 9d<p 2 )}, (3.35) 

where S(x) is given by 

{sinx if k = 1, 
X if k = 0, 

sinh x if k = — 1 . 

If k = +1, the spatial part of equation (3.35) represents a line element of a 3D hyper-surface of 4D 
ball with radius R(t) , i.e. embedded in 4D Euclidean space. If k = 0, it represents a line element of 
flat 3D Euclidean space. If k = —1, it represents a line element of a hyperbolic 3-space embedded in 
4D Minkowski space with radius R(t). 



3.3.2 The geodesic equations in the FRW universe 

To describe particles moving due to the gravitational field of the cosmological fluid, as we did in 
last section, we try to find the trajectory of the freely-falling particle, i.e. geodesic equations of our 
spacetime. We use the metric (3.35), and by applying the geodesic equation (2.7), we get the following 
conditions, i.e. the geodesic equations, which hold along any geodesic in that spacetime: 

4> = constant, 9 = constant, 

R 2 X 2 

1 H ^ — f° r massive particles. (3.36) 



R 2 X = constant, t = { R 2^2 C 



for photons. 



3.3.3 The cosmological redshift z 

In the FRW model, apart from any peculiar velocity, 20 all the galaxies are comoving, i.e. have fixed 
spatial coordinate values (x, 9, <f>). Therefore, the redshift in this model will be due to the expansion 
of the space between a stationary emitter (E), and a stationary receiver (R). Therefore, the redshift 
in this model can be written as 

vr Po(E) R{t E ) 

where po(R), and po(E) are the time component of the received and emitted photons' 4- momentum, 
respectively, and R(tu), and R(tE) are the scale factors of the universe at the reception and emission 
global time, respectively. 



3.3.4 The dynamics of the FRW spacetime 

The Einstein field equations are expressed in tensorial form as 

R^v - ^Rg^u = —kT^ - Ag^ u . (3.38) 

We assume that, on a large scale, we can model all the forms of energy in the universe, i.e. the 
cosmological fluid, as a perfect fluid. We can then make use of the energy-momentum tensor for a 
perfect fluid given by equation (2.11). We calculate the afflne connection from equation (2.4) for 
the spacetime that follows the line element given by equation (3.35), and calculate the Ricci tensor 

20 Peculiar velocity in cosmology usually refers to any deviation of the galaxies velocities from that calculated from 
Hubble's law. 
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components from equation (2.3). As a result, we get 2 independent equations for the scale factor, 21 
namely, 



R=-^(p+^R+ 1 -A c 2 R, (3.39) 
8ttG „ 1 



R 2 = pR 2 + -Ac 2 R 2 - c 2 k. (3.40) 

3 3 

Here, p stands for the total energy density of matter (m) and radiation (r), so it can be written as 
p = p m + p r . The two equations (3.39) and (3.40) are dubbed the Friedmann-Lemaitre equations, and 
if A = 0, they are dubbed the Friedmann equations. 

8ttG 

If we define the dark energy density p. such that A = — ^— p., where A is the cosmological constant, 

& 

and p total as p totaJ = p A + p = p A + p r + p m , then we can rewrite equations (3.39) and (3.40) as 

R 2 = ^P^R 2 ~c 2 k^ (p A + Pr + P JR 2 c 2 k. (3.42) 

It is common to define the normalized scale parameter a(t) as a(t) = where R(t) is the scale 

factor at any time t and Rq is its value at certain time to, which is usually taken to be the present 
time. Therefore, a = 1. 

In terms of the normalized scale factor, we can write the energy densities as 



R= — [p+^-2p A )R, (3.41) 



PaW = Pa, > P m (t) = Pr(t) = J^y> ( 3 - 43 ) 



where p A0 ,p m0 and p r refer to the present time values. The Hubble parameter is also defined as 
Moreover, we can define the dimensionless densities as 

RirC r 2 k 

where i stands for A, m, or r, and f2fe(t) is called the curvature density. 
Therefore, we can rewrite equation (3.42) as 

n A (t) + n m (t) + n r (t) + n k (t) = 1, (3.46) 

which is valid at any time. Also using equations (3.43) and (3.45), we can rewrite equation (3.42) 
again as 

H 2 (t)=H 2 {il A _ + ^ + ^ + ^fy (3.47) 
where H , £l A , Q m , f2 r and Q k stand for the values at the present time. 

3.3.5 Friedmann Models 

The analytical solutions for equation (3.47) with a zero-cosmological constant, i.e. ft A — 0, are called 
the Friedmann models. Here, we present the solutions of the dust-dominated universes, i.e. Q r =0. 
In this case, equation (3.46) is reduced to 

n ro (t)+fi fc (t) = l. (3.48) 

Therefore, using equations (3.47) and (3.48), we can write 

a 2 (t) = a 2 {t)H 2 {t) = Hi ^ + 1 - n m>0 ) . (3.49) 

21 For more detail about the derivation the reader is referred to [6]. 
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This equation can be written in integral form as 

1 f a 



' II ./, U. • M <> ),• ' 



1/2 



(3.50) 

Here, we can separate three possible cases for the universe's geometry: flat, closed and open universes. 




k = -1, open univese 



k = 0, flat univese 



= +1, closed universe 



big bang 



big crunch 



Figure 3.5: The evolution of the scale factor in dust-dominated universes 



• Flat dust- dominated universe. 



For such a case, we have k = 0, thus fl k = and from equation (3.47), Sl m D = 1, so equation 
(3.50) is reduced to 

1 f a /3 \ 2 / 3 

t = — x^dx => a(t) = i^-H B tj . (3.51) 

This case is commonly called the Einstein-de-Sitter (EdS) model. 
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• Closed dust- dominated universe. 

For such a case, we have k = +1, thus fl k < and from equation (3.47), f2 m > 1, so we solve 
the integral of equation (3.50) by taking the substitution x — ^ lm '°_ 1 sin 2 (?/;/2), where ip g [0,7r] 
and is called the development angle. This gives 

' = 2H (oj- 1)3/2 ^ " sin ^ «(*) = 2(^;°-l) (1 - C0S ^ (3 ' 52) 

This shows that the graph of a(t) with i is cycloid, i.e. a(t) starts to increase until a maximum 
value and then starts to decrease. 

• Open dust- dominated universe. 

For such a case, we have k = —1, thus fl k > and from equation (3.47), Sl m D < 1, so we solve 
the integral of equation (3.50) by taking the substitution x = sinh 2 (^/2), we get 

1 = 2 g (^";-i)3/ 2 (sinh ^ " ^ a(<) = 2(n"7°-i) (cosh ^ - 1} ' (3 - 53) 

Figure 3.5 shows, schematically, the evolution of the scale factor a(t) in the three cases discussed 
above. 
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Chapter 4 

A Discrete Model for the Universe 



4.1 Introduction 

In this chapter a model for our universe is presented. In this model the matter content is assumed to 
be discrete; identical 1 spherically symmetric islands uniformly distributed in a regular lattice. This 
attempt was first introduced in 1957 by Lindquist and Wheeler (LW) in a seminal paper [9], where 
the dynamics of a closed dust-dominated universe were studied. This model was then extended to 
investigate the optical properties for the fiat case of such a universe in another wonderful paper by 
T. Clifton and P. Ferreira in 2009 [2], The recent investigation showed a deviation in the optical 
properties for the flat case of the universe from the corresponding case in the FRW cosmology, i.e. EdS 
model. 

We start by presenting the LW model and how the dynamics for such a set-up were studied in a 
closed universe, and then proceed by introducing how this model was generalized in [2] to investigate 
the optical properties, and why this generalization is only suitable for the flat case. 

4.2 Lindquist Wheeler model for the universe 

This was an attempt to get the dynamics of a dust-dominated closed universe in the discrete case. In 
this section we briefly cover the main points of this model; the reader is referred to [9] for more details. 

4.2.1 The model construction 

From the FRW cosmology, as discussed below equation (3.35), the spatial part of a closed universe's 
line element represents the geometry of a 3D hyper-surface of a 4D ball (the comparison hyper-sphere) 
embedded in a 4D Euclidean space. Similarly, in the LW model, we take the universe's matter content 
to be discrete islands symmetrically distributed on the 3D hyper-surface of that 4D ball. As shown 
in [3], Table 1(h), it is only possible to symmetrically place N points on that surface for some values 
of N, namely, N = 5,8, 16,24, 120, or 600. Therefore, in the LW model, our universe will simply be 
N discrete islands (cells), each centered at one of these N points (the lattice points) with point mass 
m = M/N, where M is the total mass for the universe. In other words, the matter content for each 
cell will be modeled as a point mass placed at the center of that cell and everywhere else in the cell is 
empty. 

On the 3D surface, a test particle is assumed to belong to an island centered at point P if its 
closest lattice point is P. Between each pair of adjacent cells, there is a 2D array of points that are 
equally far from the centers of both cells. This array defines the 2D boundary or interface between 
the two cells. For each particular value of N, the lattice cell takes the shape of a certain curved shape 
of a regular polygon as show in Table 4.1. Inspired by the success of the approximation methods 
introduced by Wigner and Seitz in solid state physics, namely, approximating regular polygons in flat 
space by spheres of the same volume, 2 Lindquist and Wheeler suggested approximating these curved 
polygons, i.e. the lattice cells, by 3D balls. It's shown in [9] that geometrically approximating the 
curved polygons by a sphere is more accurate than approximating the flat space cube by a sphere as 

We take this case for simplicity. 
2 For more details about this approximation method, the reader is referred to [10] Chapter 9. 
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assumed and showed great success in solid state physics. Otherwise stated, a curved regular polygon 
has closer geometric properties to a sphere of similar volume than the corresponding regular polygon 
in the fiat space of the same volume. 



N 


Name for the corresponding flat 
space polygon 


5 


tetrahedron 


8 


cube 


16 


tetrahedron 


24 


octahedron 


120 


dodecahedron 


600 


tetrahedron 



Table 4.1: Closed universe cell shapes 



The Schwarzschild solution, as we discussed, gives a line element (3.3) for the spacetime in the 
vicinity of a uniform distribution of material, that line element shows that the geometry of this space- 
time is spherically symmetric but not uniform. The non-uniformity can be seen easily from the fact 
that for such line element, the Reimannian tensor Rijki is not the same everywhere. However the Ricci 
tensor Rn- vanishes everywhere except at the point mass where it is infinite. 

In [!)], Lindquist and Wheeler replaced the spherical uniform cells on the surface of the comparison 
hyper-sphere by Schwarzschild cells. 3 This means that we remove the uniformly curved spherical cells 
and replace them by non-uniform but still spherically symmetric cells. 

Because of the non-zero normal derivative 4 of the Schwarzschild potential on the cell boundaries, 
the masses from the two sides of the boundary will accelerate, in general, differently to nullify the 
difference in the normal derivative of the potential on the two sides. This leads to a relative motion 
between the lattice boundary and the lattice central mass, which is the main reason for this model to 
have dynamics. 



4.2.2 The mathematical description for the model 

We choose the Schwarzschild cells with total volume equal exactly to the total volume for the 3D 
surface of a 4D ball, i.e. 2tt 2 R 3 , where R is the radius for that 4D ball. This can be mathematically 
written as 

N 2ir ' [ ' 

where ip N is the angular separation between the two lines connecting the center of a lattice cell and 
its boundary to the center of the 4D ball. 

The radius of the lattice cell, a, and the radius of the 4D ball, R, are related by 

sin^ N = |. (4.2) 

The other necessary condition is that the new lattice cells (the patches) should be tangent to the 
surface of the comparison hyper-sphere. We can formulate that as follows; if we imagine making a 
measurement for a great circle circumference on a lattice cell, i.e. we make that measurement on a 
great circle on the cell boundary, and the circumference of infinitesimally smaller circle. Then the 
following condition should be satisfied: 

1 d(circumference) , , , „. 

n — „ , V=cosVv, (4.3) 

2n rf(radial distance) 

where ip N is as defined in equation (4.1). 

3 A Schwarzschild cell is simply a sphere but with non-uniform curvature. In other words, the geometry for such cell 
is the Schwarzschild geometry for a point mass m centered at that lattice cell center and of the same radius as that 
spherical cell it replaces. 

4 In general, the potential's normal derivative measures the rate with which a test particle moves toward/away form 
the cell boundary. 
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N 

Figure 4.1: The dependence of The ratio between the maximum radii for the comparison hyper-sphere 
in Lindquist-Wheeler universe, R , and Friedmann universe, R , on the number of the cells N for 
the closed universe. 



To find the dynamics in the model, Lindquist and Wheeler replaced Schwarzschild time coordinate, 
t, and Schwarzschild radial distance, r, in the vicinity of the boundary by a new set of coordinates 
r(r,t) and p(r,t). They constructed p(r,t) such that the boundary at any time is defined by p(r,t) = 
constant, and then they constructed r to be perpendicular to p. By formulating the conditions (4.1) 
and (4.3) in terms of the new coordinates, in [9] it was proven that the dynamics of the lattice cells 
can be described in terms of the comparison hyper-sphere radius, R, as 



^ = ^n ± 

R 2 #5sm 3 ^ R2 ' 
Therefore, the maximum radius for such universe is given by 

R LW = , . 4.5 

max iVsin 3 ^ v ' 

Here, M — mN, where M is the total mass of the universe. These results showed that if N = 600 the 
maximum radius for the 4D ball in this model and the corresponding maximum radius in the FRW 
cosmology, namely, i?^ ax = ^jr, agree to a great extent, i.e. the relative difference is about 1.2%. The 
ratio between the maximum radii is shown in Figure 4.1 for the values of N given in Table 4.1. 

The global dynamics in both cases are given by a similar evolution functional forms, and they almost 
agree on the maximum radius of the comparison hyper-sphere, in the limit of large N. Therefore, the 
large scale dynamics of the two cases of the universe is the same for large N. 



4.3 Extension for LW model to calculate the redshift 

As discussed in the previous section, LW developed their model for the case of closed dust-dominated 
universe only. They showed that this case has similar global dynamics, in the limit of a large number of 

5 This ratio is independent of the total universe mass M. 
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discrete islands, as the solution of Friedmann equation in the case of closed dust-dominated universe, 
which we discussed in the previous chapter. In [2], the main goal was to generalize such an idea for 
all the possible cases of the universe, namely, flat, closed and open and find a suitable global time 
coordinate to investigate the optical properties for such a set-up. 

The Schwarzschild coordinates can be used to describe the interior of each cell very well, however, 
they can't be used for the whole lattice, 6 because they don't match at the cells boundaries. Otherwise 
stated, the Schwarzschild coordinates intersect at the boundaries and do not overlap as it is required 
to get a smooth 3D surface on the comparison hyper-sphere surface, and as a result, at the boundaries 
there will be points with two different time coordinates. 

As discussed in the previous section, LW developed a new time coordinate suitable to match the 
time coordinates between adjacent cells at the boundaries. Since we need to investigate the optical 
properties for such a set-up we need a time coordinate which is suitable to be used not only at the 
cell boundaries but also inside the cells. To do so, the proper time of a radially out-falling particle 
was shown, in [2], to be suitable to be used as a global time for the whole lattice surface, which also 
matches at the cell boundaries, and that can be done by using the following transformation: 

/— JE - 1 + 2m Ir , N 

dr = VEdt - ¥ . 1 dr. (4.6 

1 — 2m/r 

Here, r is the proper time along the trajectory of radially out-falling particles, r and t are the 
Schwarzschild coordinates and E is a constant > 0. If we substitute the transformation (4.6) in 
the Schwarzschild line element (3.3) , we can write 



,o 1 ( 2m \ , n 2 I 2m , , dr 2 9 9 , „ , , 

ds 2 = - 1 )dT 2 + —JE - 1 + drdr r 2 d9 2 - r 2 sin 2 Odcf. (4.7 

E \ r J E V r E 

Along the radially out-falling particle trajectory, ds 2 = dr 2 and d6 = d<f> = 0, which gives 

dr\ 2 2m 



, , E-1+ . (4.8) 

dr J r 

Therefore, the 4-velocity along such trajectory is given by 8 

«*= (i,Je-1 + — ,0.li I . (4<M 




Using the metric given by equation (4.7), we can calculate the covariant components of the 4-velocity 
given by equation (4.9) to find that Ui — (—1,0,0,0). Therefore, we can easily prove that the surfaces 
of constant r are perpendicular to the trajectory of the radially out-falling particles, by taking any 
arbitrary 4-vector on a surface of constant r, which is given by 

n* = (0,n r ,n e ,n*). (4.10) 

Which follows that 

n l Ui = 0. (4.11) 
Therefore, we can use r as a suitable global time for the whole lattice (Universe). 



4.3.1 The universe construction 

Before discussing the construction, we discuss how we glue the cells of such a universe together and 
the geometry of these cells, we should also specify which regular polygon they will correspond to; in 
the two cases of open and flat universe the possible shapes are given in Table 4.2 and for the closed 
universe case the cells are curved and the corresponding flat cell shapes are given in Table 4.1. This 
problem was discussed in details by H. Coxeter in [3]. 

To construct such a set-up, we follow the following steps: 

6 If the universe in closed the lattice is corresponding to the 3D surface of the comparison hyper-sphere. 
7 We take system of units such that c 2 = 1. 

8 We taking the positive square root because we consider the out-falling particle. 
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Universe 


N 


Name for the corresponding flat 
space polygon 


flat 


00 


cube 


open 


oo 


cube 


open 


oo 


dodecahedron 


open 


oo 


Icosahedron 



Table 4.2: Flat and Open universe cell shapes 



1. We choose one of the regular lattices with cell shapes given in Table 4.1 or 4.2. 

2. At a certain initial time we place a mass m at the center of each cell. 

3. Replace each regular polygon by a sphere, and as a consequence, the lattice will have overlap and 
no man's land regions. Figure 4.2 shows how we can do that for a 2-dimensional flat universe. 9 




4. We use the same idea used by Lindquist and Wheeler in [9], and replace the line element inside 
these spheres by the line element given by equation (4.7), which is equivalent to say that we 
replace these uniformly curved spheres by non-uniformly curved spheres, i.e. with Schwarzschild 
geometry. 

5. Now each sphere will have shells which are out-falling radially according to equation (4.8) with the 
4- velocity given by equation (4.9), so we can glue the sphere at constant r, and the orthogonality 
given in equation (4.11) will guarantee the tangency between these spheres and the space in 
which they are embedded in. However, at the overlap and the no man's land regions there will 
be no tangency but the tangency will exist on average as discussed in the previous section. 

We have shown how to construct the universe, and we now present the dynamics for these cells 
and show how, by varying the value of E, we get three cases for the universe, namely, flat, open and 
closed universes. 



4.3.2 The dynamics of the cells 

According to the construction we discussed above for the universe cells (island), each spherical cell is 
composed of a continuous series of shells, and each shell has a particular radius, r(r), at every global 
time, r. We will refer to the radius of the inner shells by r(r), and the radius of the outermost shell 
(the boundary sphere) by a(r), i.e. the maximum value that r can take is a at the same global time, 
r. Therefore, we use the radius of the boundary sphere a at a particular global time r to determine 
the cell volume at this time, r. 

9 It was indicated in [0] that there will be no region common to three cells if we construct them under the condition 
(4.1). 
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Inside each cell, a particle moving with a given shell is radially out-falling according to equation 
(4.8). If we study the change in the boundary sphere radius, i.e. replace r(r) by a(r) in equation (4.8), 
we can separate the following cases. 

1. If E < 1; the solution of equation (4.8) gives that the radius of the boundary sphere increases 

2m 

until it reaches a maximum value at a max = — , and begins to decrease again, i.e. the cells 

1 — E 

start to collapse. This case corresponds to the closed dust-dominated universe solution discussed 
in Chapter 3. 

2. If E > 1; the solution of equation (4.8) gives that the radius of the boundary sphere increases 
forever and the rate of increase when r — > oo is positive, i.e. the radius will reach oo with a speed 
greater than zero. This case is corresponding to the open dust-dominated universe solution 
discussed in Chapter 3. 

3. HE — 1; the solution of equation (4.8) gives that the radius of the boundary sphere increases 
forever and the rate of increase when r — > oo is zero, i.e. the radius will reach oo with speed equal 
to zero. This case is corresponding to the flat dust-dominated universe solution (EdS model) 
discussed in Chapter 3. 

The three cases have the same schematic plot, shown in Figure 3.5. 

This generalization is suitable for the flat dust-dominated universe only. Two main reasons for 
saying that; first, as discussed in [9], the scale for the closed dust-dominated universe 10 depends on the 
number of cells N, so if we take the flat case, the spatial part will be the Euclidean 3-space and the scale 
will be the same in both cases, i.e. the radius reaches oo with speed zero. Second, and most important, 
our coordinates (t, r, 8, (f>) are suitable for the Euclidean space only, that is because, for instance, in 
closed case, they will not cover all the space, i.e. in the vicinity of the maximum expansion, they will 
not cover all the no man's land regions. Therefore, from now on we shall take E = 1. The solution 
for equation (4.8) in the flat universe case, assuming 2m = 1 for simplicity, is r — r = |(r 3 / 2 — r^ 2 ). 
Figure 4.3 shows how the radii for different shells and for the boundary sphere increase with r. 




radial distance 

Figure 4.3: Shows the solution for equation (4.8) in the flat case, i.e. E = 1. Here, we take 2m = 1 and 
assume that the radius of the boundary sphere is a — 2.0 (the purple curve) at a certain initial time 
t =0. At t — 0,The inner shell will be smaller initial radii, and in the figure we show the change in 
different shells' radii with r 



'We mean here by the scale the maximum radius of the universe. 
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4.3.3 The geodesic equations in terms of the new coordinates 

We derive the geodesic equation using the Lagrangian method, as discussed in Chapter 3. Therefore, 
using the line element given by equation (4.7) with E = 1, our Lagrangian can be written as 

L = g^iTiP = - ^) f 2 + - r 2 - r 2 9 2 - r 2 sin 2 9 ft. (4.12) 

Here, the dot refers to the derivative with respect to the parameterization parameter A used to param- 
eterize this geodesic. Using the Euler-Lagrange equation — [ — — ) = — — , we can write that along 

dX V ox^ I ox^ 



the geodesies the following relations always hold: 

(t(1 - 2m/r) + ry/2m/r) = 0, (4.13) 



d 

dX 



^-^^fr = - m ,T 2 +r9 2 +r sin 2 9 ft, (4.14) 
d\ dX r 1 

J- (r 2 e^ = r 2 sin 9 cos 9 ft, (4.15) 



d (r 2 sin 2 = 0. (4.16) 



Equations (4.13) and (4.16) give 



r/A 



f(l - 2m It) + r^2m/r = B, (4.17) 



r 2 sin 2 9 ' 



(4.18) 



where B and are constants. Equation (4.15) can be rewritten as 



r 



J? cos ( 



+ 2-9-^ — s-=0. (4.19) 
r r" 4 sin 9 

d ( J 2 \ 

If we multiply equation (4.19) by r 9, we can write — ( r A 9 2 + = 0. It follows that 



d\ 



sm 



sin' 9 

where J is a constant. Therefore, we can write that 

J 2 Jl 



r*9 2 + ^l- = J 2 , (4.20) 



Q2 



r 4 r 4 sin 2 i 



(4.21) 



The 4-velocity for a particle moving geodesically is determined by the functions f , r, 0, and <fi which 
obey equations (4.13) - (4.16), equations (4.21) and (4.18) determine 9 and <j), respectively. Now, as 
we did before, we can separate the two case of photons and massive particles. 

1. Photons: To get the full photons trajectory, instead of using equation (4.14) and (4.17) to 
determine f and r, we use the null-constraint g^i^i^ = with equation (4.17). We can easily 
show that 

Therefore, equations (4.17), (4.18), (4.21) and (4.22) determine the geodesic of a photon in the 
cell. 
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2. Massive particles: Similarly we use the null-constraint g^x^x^ = 1 with equation (4.17), we can 
easily show that 

9 2m , J 2 ( 2m\ 

r 2 = 1 + B 2 - 1 . 4.23 

r r A \ r J 

Therefore, equations (4.17), (4.18), (4.21) and (4.23) determine the geodesic of any massive 
particle in the cell. 

We can also prove that the trajectory of massive particles moving with the bounding spheres is 
a geodesic. The geodesic of a massive particle with J = = 0, i.e. 6 = <j> = and B = 1, i.e. 
r = 1, so that f = yjlrajr. Which has the same 4- velocity, given in equation (4.9), of a massive 
particle that is radially out-falling with a spherical shell. Therefore, such a particle is moving 
geodesically. 11 

It is also worth noting that we did not restrict ourselves to the equatorial plane as we did in Chapter 
3, because we will need to match the coordinates of adjacent cells at the boundaries, and so rotating 
the coordinates would make the matching process too complicated. 

Now, we have a full description for the photons trajectories in each cell. To investigate the optical 
properties for the universe, we need to know how the photons' trajectories behave when the photons 
move from one cell to another. One might think that this can be easily done by making a transformation 
for the spatial coordinates between the two cells. However, in [2], it was shown that this is not true, 
because the null-constraint in the first cell, using such transformation, will give an equation in terms 
of the second cell coordinates, and this equation contradicts with the null-constraint equation of the 
second cell. 



4.3.4 Matching the coordinates at the boundaries of the cells 

The main idea here is that if an observer at a cell boundary measures the photon frequency or direction, 
at a particular global time, using the coordinates of both cells, then, in two measurements, he should 
get the same results. 

When photons pass from one cell to another, most of them will pass through the overlap or the no 
man's land regions. The photons, that will cross at intersection points, will be a set of measure zero. 
As shown in Figure 4.4 on the left part, while a photon is crossing from the boundary of cell 1 to the 
boundary of cell 2, i.e. moving between A and B, the distance A-B changes, and therefore, we expect 
a change in the frequency. The same thing happens if photons pass through the no man's land region, 
the frequencies in both cases will change in an opposite sense and on average we expect that when the 
photons cross a boundary, their frequency will not change. Otherwise stated, if the photon in cell 1 at 
the boundary has f = fx, and at the boundary of cell 2 has t = f-z, on average, we have fx = T2- This 
approximation has been shown to be valid numerically in [2]. 




Figure 4.4: On the left, we see two adjacent cells and two observers, A and B, on their boundary. On 
the right, it shows a cell with an observer on an inner shell, C, and another one on the boundary, D. 

The 4-velocity for the photons' trajectories can be decomposed as 

k a = (T,r,e,<j)) = (-u b k b )(u a + n a ), (4.24) 

where u a is given in equation (4.9) and n a is given in equation (4.10). It follows that 

k a = (f, r, 6, 0) = f(l, y/2m/r + n r , n e , n*). (4.25) 

11 This is another reason for making this calculation for the flat universe; the particle moving radially with the shells 
are not moving geodesically in the other cases of the universe. 
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This means that in any cell, the photons' trajectories read 



rrr , 

: fn*, (4.26) 
rn + Ty/2m/r. 



Therefore, to transform a photon trajectory from one cell to another, we follow the following steps: 

1. At the boundary of cell 1, we take t x , r 1 , 9\ and and use equation (4.26) to determine n r . n® 
and nf. 



2. Determine n a for cell 2 from n a — Tr4-n b , where x a and x a are related by 12 

2 2 (9t 1 12 J 



K = r t + X o ~ 2r i cos sin 9, 

9 (r. cos(i), sin 6> — Xr,) 2 

cos 6 — — — - — (4 27) 

2 A 2 + r 2 sin 2 X - 2r 1 A cos 1 sin 9 ± ' 



cos 2 1 



r 2 cos 2 i 



Xl + r\ sin 2 9 X - 2r x A cos 4> 1 sin ( 



where x a , or x b , refers to r, 9 and 

3. Given that at the boundary on average t 1 ~ f 2 , we use the results of step (2) and equation (4.26) 
to calculate r 2 , 9 2 and <p 2 . 

4. We use these values as initial values to propagate the photons in cell 2, and upon reaching the 
boundary of cell 2, we take r 2 , r 2 , 9 2 and <f> 2 and repeat the above steps. 

It is also worth noting that inside each cell the shells are out-falling according to equation (4.8), which 
shows that for the two observers C and D moving geodesically. 13 For the observers shown in the 
right part of Figure 4.4, there will be relative motion between them because they have different radial 
coordinate, and therefore different radial velocity, i.e. there is a relative motion between C and D. This 
relative motion will lead to a redshift inside each cell. 

As discussed above, on average, there will be no changes in the photons' frequencies at the bound- 
aries, which means that the redshift, in this model, occurs only inside the cells and comes only from 
the fact that the shells that are radially out-falling follow the 4- velocity given in equation (4.8). Now 
that we have the fall trajectory of the photons in the universe, i.e. trajectories along which the photons 
travel between any two point in two different cells and across any number of cells, we can calculate 
the redshift occurs in this universe. 



4.3.5 The redshift 

We use the observers moving geodesically with the radially out-falling shells in the cell as the closest 
analog to the co-moving observers in an EdS model to compare the redshift in the two cases, namely, 
the flat matter-dominated discrete universe and the perfect fluid universe (EdS model). 

The photon's frequency, emitted from, or received by, such observers, is given by f = v = —u a k a , 
where u a and k a are given in equations (4.9) and (4.25), respectively. Therefore the redshift can be 
expressed, in general, as 

v e f, , 
l + = (4.28) 

where r e and r r are the values of f along the null-trajectory when the photon was emitted and 
received, respectively, and similarly for v values. In general, tracking how f changes inside all the 
cells, the photon crossed, is only possible numerically. However, we present here an approximation to 
calculate it analytically in terms of the corresponding redshift in the perfect fluid universe. 

12 Using the transformation x 2 = x x — Xo, y 2 = y 1 and z 2 = z x , we find the spherical polar coordinates in cell 2 in 
terms of the spherical polar coordinates of cell 1. Here, Xq is the cell width at the global time when the photon crosses 
the boundary. 

13 The reason for saying that they moving geodesically is discussed below equation (4.23). 
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* An analytic approximation to calculate the redshift 

First, we try to find the redshift inside a single cell k, 1 + 6Zk- 
Equation (4.22), using equations (4.20) and (4.18), can be rewritten as 

r 2 = B 2 - 1 ~^™ /r (r i 9 2 + r 4 sin 2 9<j) 2 ) = B 2 - (1 - 2m/r)r 2 (6 2 + sin 2 9<j> 2 ). (4.29) 

Assuming that r 3> 2m, we can write that B = y r 2 + r 2 (9 2 + sin 2 9<f> 2 ). Therefore, in this approxima- 
tion, B expresses the total velocity of the photons inside each cell. We will assume that it is constant 
for all cells. 14 

Moreover, from equation (4.17) we can write that f = B 1 ; where a = -g, the fraction 

of the radial velocity (r = v = y/2m/r) with respect to the total velocity. If we take the limit that 
r ^> 2m, i.e. v is small, we can rewrite f as 

f = B 1 ^™ ~ B{1 - av)(l + v 2 ) ~ 5(1 - u(a - v)). (4.30) 

1 — tr 

Now we can calculate the redshift of a given photon inside a single Schwarzschild cell. We take f in 
and r out to be the values of f when this photon enters and leaves the cell, respectively. And since 
the photon enters and leaves any spherical cell at the same radius, 15 v in ~ v out ~ v ~ \j2mj a{r). 
Therefore, the redshift of a photon while it is crossing the cell k can be expressed as 

• 1 — v k (a k n — v k ) 



1 + 5z k = 



j-k I _ v k( a k _ y k\ 

out V out ' 

[1 - v k (a k n — v k )][l + v k (a k M — v k )} ~ 1 + v k (a k ut — a* 



(4.31) 



Second, we calculate the total redshift. To do so, we imagine that we have three points 1, 2 and 3 on a 
photons' trajectory, and that we place observers to measure photons' frequencies, V\, and V3 at these 
points. The redshift between 1 and 2 is 1 + zi_2 = The redshift between 2 and 3 is 1 + Z2-3 = 
The redshift between 1 and 3 is 1 + Z1-2 = ^ = ^Sf = (l + z i-2)(l + ^2-3)- Therefore, for a photon 
crossing n-cells, the total redshift is 16 



1 + z = Y[(l + Sz k )~l + ^2Sz k . (4.32) 



fc=i fc=i 



Here, we will make a further assumption, namely, we assume that r will not change much inside each 
cell, and therefore, r k = —r k =>■ a = — a k =>■ d fe — d* 1 = 2a k . This means that, using equation 

' ' m out m out out in ' *-* 1 

(4.31), 5z k = 2a k v k = 2 v / 2m/a k a k . Substituting this in equation (4.32), we can write that 





1 + ^ = 1 + 2^/— a k ^l + 2j— / —^a k . (4.33) 

k=l 

Here, as discussed in Chapter 3, we used the fact that for a flat dust-dominated universe air) cx r 2 / 3 . 
As we assumed before, we ignore the expansion of the cells while a photon crosses them, so we can 
write = — 2a k . This gives, ^ ~ — 2a k — —2a k (r /t ) 2 / 3 . If the photon crosses cell k non-radially, 
i.e. only crosses j3 k from the cell width (2a), we write that ^ ~ —2a k /3 k (T/T ) 2 ^ 3 . Substituting this in 
equation (4.33) gives 



1-!-: : l-,/^-n-r)- = l-^T /^-\ln /r ' 



I 2m T o r ( a 

a o a J \/3 k 



Ik 



Tr 



(4.34) 



a. 



(7) In ( J ) = 1 + (7) ln(l + z FRW ) ~ (1 + z F 1 



a k 



Here, (7) = ^ ^ is the ensemble average. 

In [2] (Appendix E), the value of (7) was approximated geometrically to be 2/3, however, numeri- 
cally, its value was estimated, also in [2], to be 7/10. 

14 This is consistent with the constancy in the light speed in the perfect fluid universe. 

15 Here, we assume that the radius of the cell will not grow much while the photon is crossing the cell. 

16 We are assuming that the redshift inside the cells is small. 
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4.4 Discussion 



The value of (7) was estimated both analytically and numerically in [2] and was shown to be approx- 
imately of the same value. This showed that there is a deviation in the value of the redshift in the 
discrete case from the corresponding values in the perfect fluid case. This means that the averag- 
ing process of the spacetime, i.e. taking perfect fluid assumptions does not in general commute with 
the evolution of the universe. Moreover, this deviation shows that we should take the effect of the 
inhomogeneity in our interpretation of redshift and luminosity distances. 

In the Lindquist and Wheeler model and in the extension trail we discussed above, only the case 
of equal masses regularly distributed was considered. However, the current surveys show that the 
mass distribution in the universe is due to a stochastic process during the formation of the universe 
structure, which means that considering such a regular distribution is a huge idealization of the universe 
and therefore, the next step to improve this model would be taking irregular mass distributions and 
considering islands with different masses. This would make the calculation very difficult but as we 
referred above the averaging process does not commute with the universe evolution, so we expect a 
further deviations from the regular case considered here. 

It is also worth noting that this model depends mainly on the assumption of approximating the 
spacetime geometry inside each discrete island by a Schwarzschild geometry which can be seen as a 
weak point in this model. 
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